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Abstract 

We show that the exact ground state for a class of extended Hubbard models 
including bond-charge, exchange, and pair- hopping terms, is the Yang "77- 
paired" state for any non-vanishing value of the pair-hopping amplitude, at 
least when the on-site Coulomb interaction is attractive enough and the re- 
maining physical parameters satisfy a single constraint. The ground state is 
thus rigorously superconducting. Our result holds on a bipartite lattice in 

any dimension, at any band filling, and for arbitrary electron hopping. 
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Interest in itinerant, strongly interacting electron systems has exploded since the dis- 
covery of high-T c superconductors, for which the interplay between itinerant magnetism 
and insulating behavior is believed to play a crucial role. The Hubbard model |1] provides 
the simplest description of such systems, by assuming that the itinerant electrons interact 
only via an on-site term, and it has therefore been intensively studied. The original Bethe 
Ansatz solution of the one-dimensional model at half-filling has in recent years been 
supplemented by a number of additional rigorous results (see |§ for a current review). For 
instance, for the attractive (i.e., "negative U") Hubbard model, long studied as a model 
believed to have a superconducting ground state recent articles have established the 
existence of "off-diagonal long-range order" (ODLRO 0), but these results have been re- 
stricted either to U — > — oo || or to bipartite lattices in which the number of sites on one 
sublattice is not equal to that on the other Indeed, although it is known that 

the Hubbard Hamiltonian has certain eigenstates (the "^-paired" states) that have non-zero 
pairing and ODLRO and hence are superconducting, Yang has shown that these states can 
never be the ground state for the pure Hubbard model on a standard bipartite lattice ||. 

Apart from the Hubbard model itself, "extended" Hubbard models have also attracted 
considerable recent interest, in part because it has been recognized PH|JTT| ] that the addi- 
tional interaction terms (discussed but eventually neglected in the original papers [|l||) could 
be relevant in stabilizing some novel (e.g., ferromagnetic or superconducting) phases. For 
instance, Strack and Vollhardt Jl2| proved rigorously that the ground state for a class of 



extended Hubbard models is a saturated ferromagnet at half-filling for any non-vanishing 
value of the exchange term, provided that the Coulomb repulsion is strong enough. 

In this Letter we show that in appropriate regions of the parameter space of a class of 
extended Hubbard models, the superconducting phase is the stable ground state for any 
non- vanishing value of the pair-hopping term fl0| , [T3|j . More precisely, by simultaneous use 
of Yang's states and of a suitable generalization of Strack and Vollhardt's techniques, we 
prove rigorously that the ground-state of the extended Hubbard model at any band filling 
and for average magnetization m = is superconducting, provided that the on-site Coulomb 



interaction is attractive enough and that the remaining physical parameters satisfy a single 
reasonable constraint. Moreover, by a particle- hole transformation, we map our result into 
a sufficient criterion for stability of ferromagnetism in the repulsive case at half-filling and 
for arbitrary m, which criterion is more general than the one given in |H2| . Indeed, this 
relation between superconductivity for negative U and magnetism for positive U can readily 
be understood in terms of the 5*0(4) symmetry [14| underlying the model. Importantly, in 



contrast to other recent work ||15|| , the supersymmetric condition t = X (see (pj below) is in 
general neither required nor fulfilled by our superconducting ground state. 
Our extended Hubbard Hamiltonian reads [HI: 



H = ~t Z) £ C L Ck >° + U Y, ^j.T^ja + X Z Y,( n i,-° + ™k,-a)c}, CT C k,<x 
<j,k> cr j <j i k> a 

+ -o J2 n i n * + Y J2 J2 c L c U' c j,<x' c k, CT + -k J2 Y, c lo c l-° c *-° c ^° > C 1 ) 

<j,k> <j,k> a,a' <j,k> <? 

where c^ a ,Cj >a are fermionic creation and annihilation operators ({cj )(T ', Ck )Cr } = 0, 
{ c j,o-' c k,o-'} = ^j,k ^cr, o-'I, = c j,o- c j,o-' n j = Scr n j,o- ) on a d-dimensional lattice A (j, k G A, 
o G {"IT, JJ-})) an d < j, k > stands for nearest neighbors (n.n.) in A. In (pj) the first term 
represents the band energy of the electrons, and the remaining terms describe their Coulomb 
interaction energy in a narrow band approximation 0: U parametrizes the on-site diagonal 
interaction, V the neighboring site charge interaction, X the bond-charge interaction, W 
the exchange term, and Y the pair-hopping term. An explicit evaluation of the relative size 
of these contributions - all generated from onsite and n.n. matrix elements of the Coulomb 
interaction - was already given in [pj. It is worth emphasizing that most of the following 
analysis can be extended in a straightforward way to the case in which the interactions in 
(|l|) are not confined to neighboring sites. However, to avoid cumbersome notation, we have 
chosen to limit ourselves to the present case. 

The exchange and pair-hopping terms in (fl|) can be written more conveniently in terms 
of the conventional spin and pseudo-spin operators and S^ a \ a = x,y, z, 



^ +) = H'Mtii ' ^ = (t^ . *J W = ^Kt + n JlX - 1) (2) 
sf } = sf } ± iSi y) , sf^^iis^ . 

which are known to generate two orthogonal su(2) algebras. For a bipartite lattice with an 
even number of sites, the second line of ([I]) reads then 

2V £ s^s^-w £ s r s k -y £ Sj-s k + c , (3) 

<j,k> <j,k> <j,k> 

1 l.W 
where V — V — ~ Y) an d C = ~(V —)qN(2n — 1), with g number of nearest 

neighbors in A, N number of sites, and n average electron number per site. 

H can be easily seen to commute with su(2)f, generated by S a = y~] S?, a = x,y,z. 

j 

This observation is in fact relevant in recognizing the (N + l)-fold degenerate saturated 
ferromagnetic state \ip m >, 

(n+m)iV 

\iff m >= (5+) 2 \H . . .|> , |m| < min{l,n} , (4) 

as eigenstate of if at n = 1 (half-filling) and average magnetization per site m. Analogously, 

one can construct su(2) sc , generated by S a — The latter does not commute with H. 

j 

Nevertheless, it is easily checked that the states - known as 77-pairs || - defined by 

\r) n >= (£+)^|0> , 0<n<2 , (5) 

with |0 > being the electron vacuum, are eigenstates of H at band filling n with m — 0, 
provided that V' = 0. Interestingly enough, these states are precisely those which Yang has 
proved to exhibit ODLRO, which property has been shown to imply both Meissner effect 
and flux quantization [ 16| , i.e. superconductivity. Notice that the 77-pairs (||) differ from the 



pairs defined in Jl5| by a factor of (—)'•"; this will be essential in the following. 



Assuming henceforth V = 0, we want to investigate under which circumstances \r/ n > is 
indeed the ground state. We proceed by rewriting, by means of suitable operator identities, 
the hopping and bond-charge repulsion terms as sums of positive definite operators having 
zero eigenvalue on (|5]) plus contributions which simply renormalize the other terms in the 
Hamiltonian, which now reads 



h = -u'Us^] 2 -W E Sj-Sk-y E (sj • s k - i) + J E s\ z) s^ + c (e) 

j <j,k> <j,k> V J <j,k> 

E \\t - ^il7|o] k o jk + e [l* - ^ll 1 - tI^W^w + \x\{q- 3 ±,M^ + R k« R ^)] } > 

<j,k> I cr J 



where U' = 2{U + q[\t - X|(| 7 |a 2 + |1 - 7 |/3 2 ) + 4|X|]}, W = W - (3 2 \t - X\\l - j\, 

Y' = Y-\t-X\[^ + 11^1], J = \t - X|[2| 7 |(a 2 + 1) - |1 - tI(/5 2 + ^)], and C" 
1 g 

C — -(Un + Y-)N. Here a 7^ 0, /3 7^ 0, and 7 are free paramaters, and 
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Ojk = a(£« - S£>) + - (i T c jjT + cj a c ka 



P . _ 1 

^jk,a 2 



Qjk,a = — ?= (nj_ ff Cj iCr + 6»n k _ CT c k;CT ) , 

^k, ff = -^[(l-nj,- ff )qj )<T + e(l-n k ,_ ff )c k)ff ] , (7) 

with e = sg?7,[(t — X) 7 ], 77 = sgn[(t — X)(l — 7)], and 9 = sgn[X}. 

For U' < and y > the lower bound of the on-site and pair-hopping terms in is 
zero, which coincides with their eigenvalue on the states \rj n >. Moreover, 

Ojk\rj n >= Pjk,a\Vn >= Q\v,aW >= ^jk,<r|^n >= , (8) 

and - H being a positive definite form in these operators - \r] n > is the ground state also 
for the second line of (^). The freedom in the choice of a, (3, and 7 permits two of them to 
be fixed so that W = = J. In this case, one finds that the state \r] n > is the ground state 
with energy E gs = — (U — qY)nN whenever U' < 0, and Y' > 0. In fact, we can obtain an 
even larger region of values of U for which \r] n > is the ground state if, instead of fixing a, 
(3, and 7 as above, we first express both the Ising-like term and the ferromagnetic exchange 
term in (|6]) through the following operator identities, 



w ' E s r s k = \W'\ ( E ^Ik^jk - 2gE^f] 2 ) 

<j,k> \<j,k> j / 



w E s ( j z) sl z) 



<j,k> 



j, E S ( Z ) S (*) = ! A 1 E [5 w + sgn{f)s ^f - q ^[S\ z) f ] , (9) 
<j,k> V <j. k > j 



with 
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Bjk = 72 W' T ° ja " s 9 n ( W ') c U c ^l) > ( 10 ) 
and J' = J — W . The identities ([5]) - when inserted in (||) - further renormalize the 
coefficient of the Coulomb interaction term {i.e. X^j/Sj^] 2 ), so that the inequalities which 
have to be satisfied in order that \rj n > be the ground state now read 

U' + q(2\W'\ + \J'\) < , Y' >0 , (ff) 

and are still functions of a, /3, and 7. Eliminating a using the second of Eqns. QTT]) , one is 

left with a single inequality for u = — , 

qt 



u < — < 4x + |1 — x\ 



27 2 (1 - 7) 2 



iw (1-7) S 



x| y — B 



B 2 1-i 



(12) 

2/ |!-7l 



where x, y, and u; are X, Y, and W expressed in units of t, and B = 
< B < y. Equation QT2] ) can be optimized by fixing the remaining parameters /? (through 
B) and 7 variationally. The result gives a sufficient condition for the superconducting 
state \rf n > to be stable ground state. This ground state is unique for Y' > 0, for in 
that case it is the unique ground state of the pair-hopping term in The constraint 
V — (i.e. V = —(W — Y)) must be satisfied. The above rigorous result holds for any 
bipartite lattice, in any dimension, and importantly for arbitrary values of x. In particular, 
for v = x = y = w = one finds directly from ( |i~2"D that the pure Hubbard model has a 
superconducting ground state with ODLRO at least for u = —00, in agreement with 0. 

The final explicit form of the general result (|T2"D is too long to be written in the present 
letter, and will be reported elsewhere |T7| . Here we simply plot in Fig. [I] the actual boundary 



of our rigorous superconducting region in the u vs y plane at different x values. The 
inequality ([T2| ) and Fig. [I] show that any y 7^ there is a region of u values for which the 
system is superconducting, and its size increases with increasing y, at least for y < |1 — x\. 
It is quite natural that a non-vanishing value of y can stabilize the superconducting phase, 
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in that it removes the degeneracy of the expectation value of the Hamiltonian on states 
|$ n (0) >= e*^(cj |cj jj] n ~, which, for arbitrary 7^ vr and x 7^ 1, are not eigenstates of 
H. As pointed out by Yang in ||, this observation implies that \rj n >= |$ n ( 7r ) > can not 
possibly be the ground state at y — 0. On the contrary, for 2 = 1 and y = 0, |$ n (^) > is an 
eigenstate of H for all <f> and can in principle be the (degenerate) ground state. Indeed, from 



(|i2|) we see that this is the case at least for u < —4. For y 7^ 0, x 7^ 1, < $ n (<^)|i7|$ n (0) > 
becomes a function of (j), in fact minimized by <f> = it, and l^nfV) > turns out to be the 
ground state at least in the region of U values satisfying fll2l) . Finally, for x = 1 and y 7^ 0, 
there are two choices of <f> which correspond to eigenstates of H, (f> = tt (for V = -^W ~ Y)) 
and = (for V = + ^) < 0). The first corresponds to the ground state at least for 
u < — (4 + y + w) (see (|12|)), whereas the region of stability of the solution corresponding to 
the second was already discussed in [ |15|j . 

The relation (|T2|) has a solution only for negative values of u. On the other hand, the 
physics of high-T c materials suggests that the actual value of the on-site electron interaction 
is strongly repulsive. Even if the electron-phonon coupling reduces the effective value of the 



Hubbard interaction ||18|| , its sign is still expected to be positive. However, one should keep 
in mind two points. First, (|T2|) is a sufficient condition, and thus does not eliminate the 
possibility of having \n n > as ground state even when it is not fulfilled. Again, the easier 
case x = 1 helps clarify this point. There, an exact solution in 1-d [[0|] shows that ODLRO 
and superconductivity still survive as part of the degenerate ground state up to moderately 
positive values of u, which values are now band-filling dependent. We thus expect this 
behavior to persist even when the condition x = 1 is relaxed. In particular, it would be 
extremely interesting to work out an exact Bethe Ansatz solution of ([!]) in d = 1, at least 
at x = (for u = w = the latter would be the superconducting "t — Y" model, which is 
the particle-hole transformation of the ferromagnetic t — J model p(J). Second, and very 
importantly, knowing rigorously both the superconducting nature of the ground state and its 
explicit form in any dimension provides powerful benchmarks for the approximate methods 
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required to examine more realistic models. 

Apart from the superconducting solution, our expression (^|) for H allows us also to 
recognize a region of the parameter space characterized by ferromagnetic order. Indeed, it 
is easily seen that the state \ip m > of (|]) can be obtained from the state \r] n > of ([5]) by the 
following unitary particle-hole transformation: 

c j,f -» c j,t , 4,1 ~* H lj ' c j,i • ( 13 ) 

The same transformation maps Sj"' 1 into Sj , and the consequences for H as given by @ 
can be worked out directly. Let us call the transformed Hamiltonian H. For X = 0, H is 
still an extended Hubbard model, in which the on-site Coulomb repulsion term has opposite 
sign, W and Y have exchanged their roles (the first becoming the pair-hopping amplitude, 
and the second the exchange coupling), and the operators of ([?[)- (fL0|) have been redefined 
accordingly. Moreover, an arbitrary neighboring-site Coulomb repulsion term can be added 
to H, as now it simply renormalizes the coefficient J of the Ising-like term in (^). The 
discussion following (|6|) can be used to examine the conditions under which the saturated 
ferromagnetic state \ip m > is the ground state of H. A straightforward calculation shows 
that the result is identical (in form) to the one given in (|T^), apart from the sign of u, the 
inequality hence becoming a lower bound for positive u. Further, now it is the exchange 
coupling which cannot be zero in order to have a stable ferromagnetic phase. This result 



is in full agreement with ||12|| . In fact, our lower bound can easily be seen to coincide with 
expression (6) of |12| for 7 = 0, whereas it is lower than that if 7 is fixed variationally. Again, 
a more complete discussion of this case will be given elsewhere |17]]. Notice that in |12[ the 



freedom in the polarization of the saturated ferromagnet was not explicitly incorporated, 
which hid the power of the particle-hole transformation. 

In summary, we have shown rigorously that a large class of extended Hubbard models 
on bipartite lattices has a superconducting ground state for negative U and non-vanishing 
pair-hopping amplitude. The conditions derived here are sufficient, depend in a trivial way 
on the dimension, and do not depend at all on the band filling. Of course, this does not 
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exclude that a superconducting ground state can exist even for moderate positive values of 
U. If this is the case, we expect that the dimension and the band filling should become 
crucial, as happens for instance in the t = X case. Work is in progress along these lines. 
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FIGURES 

FIG. 1. The boundary of the superconducting phase in the — u-y plane for different x values, 
x = 0, .5, 1 (represented by continuous, dashed, and dot-dashed lines, respectively), for v = 0, 
w = y. 
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